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1 .  Introduction 


The  bulk  of  the  atmosphere,  the  upper  layers  of  the  oceans  and 
continents,  and  the  ice,  snow,  and  other  material  that  often  cover  these 
layers,  may  be  said  to  make  up  the  climate  system.  Processes  that  Influence 
the  system  but  are  not  in  turn  influenced  by  it  may  be  considered  external. 

If  a  climate  has  become  established,  a  subsequent  small  change  in  some 
external  condition  should  lead  in  due  time  to  the  establishment  of  a  new 
climate,  probably  differing  slightly  from  the  earlier  one.  For  example,  a 
small  Increase  In  solar  output,  represented  perhaps  by  a  0.2-degree  increase 
in  the  average  planetary  temperature  T«* — the  average  temperature  that  the 
earth  would  acquire  If  it  were  a  black  body  in  radiative  equilibrium  with 
the  sun — might  actually  cause  a  0.6-degree  increase  in  the  global  mean  sea- 
level  temperature  To,  because  of  the  greenhouse  or  some  other  effect.  In 
such  a  case  we  would  say  that  the  sensitivity  of  To  to  To*  Is  2.6  (degrees 
per  degree).  We  can  also  speak  of  the  sensitivity  of  one  quantity  to  a 
dimensionally  different  quantity — for  example,  the  sensitivity  of  globally 
averaged  rainfall  to  To*,  or  the  sensitivity  of  To  to  atmospheric  carbon 
dioxide  content. 

Under  suitable  circumstances  the  sensitivity  of  one  quantity  to  another 
may  be  infinite,  l.e.,  an  infinitesimal  change  in  some  condition  may  Induce 
a  finite  response.  Infinite  sensitivity  is  especially  likely  to  be  found  in 
idealized  models.  Fig.  1,  for  example,  is  based  on  a  low-order  moist 
hemispheric-circulation  model  which  was  developed  under  a  succession  of 
contracts  with  AFGL;  for  documentation,  see  Lorenz  (1984).  We  shall  refer 

A* 

to  the  model  as  M.  The  three  curves  show  the  values  of  the  atmospheric  sea- 
level  temperature  T  that  would  be  in  equilibrium  with  the  indicated  values 
of  the  "local  planetary  temperature"  T\  for  different  values  of  one  of  the 


controllable  parameters,  when  the  horizontal  variations  of  T*  have  been 
suppressed. 


In  the  model  M,  liquid  water  which  has  condensed  In  the  atmosphere  is 
carried  along  for  a  while  before  falling  out  as  rain;  clouds  therefore  have 
a  "half-life*  h.  For  the  left-hand  curve  In  Fig.  1,  h  has  been  sej^*tf"T8 
hours;  a  rather  high  sensitivity,  peaking  at  about  fi.O^qoedfS  within  a 


fairly  limited  range  of  T\  centered  near  267  JC.^For  the  middle  curve  h  has 
been  Increased  to  36  hours;  we  find  two  stable  equilibrium  values  of  T 
surrounding  an  unstable  equilibrium  when  T*  lies  in  a  very  narrow  range 
centered  near  276  K.  As  the  upper  or  lower  limit  of  this  range  is 
approached  from  below  or  above,  respectively,  T  becomes  Infinitely  sensitive 
to  an  increase  or  decrease  in  T*.  For  the  right-hand  curve,  where  h  is  54 
hours,  the  range  of  T*  for  which  there  are  multiple  equilibria  is 
considerably  wider. 

The  high  sensitivity  occurring  in  M  apparently  results  from  a  cloud- 
albedo  feedback  process.  In  M  the  albedo  a  has  been  made  proportional  to 
the  cloud  cover,  which  in  turn  is  parameterized  in  terms  of  the  relative 
humidity  r.  Although  no  relation  between  r  and  T  is  built  into  M,  numerical 


solutions  indicate  that  in  general  r  is  lower  when  T  is  higher  An  increase 
in  T*  will  have  the  Initial  effect  of  Increasing  T;  the  accompanying  drop  in 
r  and  hence  in  a  allows  more  solar  radiation  to  heat  the  earth,  thus 
further  Increasing  T.  With  a  high  value  of  h,  when  the  atmosphere  has 
difficulty  in  shedding  its  clouds,  the  cloud-albedo  feedback  process  is 
enhanced.  Auxiliary  computations  where  a  is  held  fixed,  although  r  is 
allowed  to  vary,  exhibit  only  finite  sensitivity  of  T  to  TV 

Fig.  2  shows  the  rate  P(T*,T)  at  which  T  will  vary  in  seeking 
equilibrium  with  a  fixed  T\  for  the  case  where  h  =  36  hours.  Actually, 


^jrfien  all  horizontal  variations  are  suppressed,  M  retains  three  dependent 
variables — the  total  dew  point  W,  the  sea-surface  temperature  S,  and  T. 

When  T*  is  fixed,  T,  W,  and  S  will  ordinarily  undergo  three  modes  of 
variation.  Two  of  these,  characterized  by  rapid  and  sometimes  oscillatory 
damping,  represent  mutual  adjustments  among  T.  W,  and  S.  The  third  mode, 
which  Is  the  one  depicted  in  Pig.  2,  is  characterized  by  slow  damping,  and 
represents  a  simultaneous  adjustment  of  T,  W,  and  S  to  T*.  In  a  typical 
time-dependent  solution  only  the  third  mode  is  noticeable  after  a  month  or 
two.  We  shall  treat  P  as  the  rate  of  diabatic  heating,  even  though  the  more 
rapidly  oscillating  modes  are  also  diabatic. 

The  zero  line  is  of  course  the  same  as  the  middle  curve  in  Fig.  1;  the 
other  isolines  are  somewhat  similar  in  shape.  The  response  can  be  very 
slow;  sometimes  many  years  are  needed  for  T  to  come  within  one  degree  of 
equilibrium.  For  example,  when  T*  -  278  K,  Just  above  the  range  where 
infinite  sensitivity  exists,  and  T  Is  initially  low,  say  250  K,  T  will  at 
first  appear  to  be  slowing  down  toward  an  equilibrium  at  about  268  K.  Only 
as  T  nears  268  K  does  it  become  apparent  that  T  is  not  about  to  settle  down 
there,  and  in  due  time  T  rises  more  rapidly,  finally  slowing  again  and 
approaching  a  true  equilibrium  at  306  K. 

There  is  no  compelling  evidence  that  the  cloud-albedo  feedback  process 
encountered  in  M  operates  in  a  similar  manner  in  the  real  atmosphere. 
Possibly  some  other  real-world  processes  operate  in  this  manner.  In  this 
report  our  interest  will  be  in  any  process  which  can  produce  a  high  or 
perhaps  infinite  sensitivity  under  a  restricted  set  of  external  or  Internal 
conditions. 

Pig.  3  shows  the  variations  of  To  for  a  400-year  numerical  run  made 
with  M,  starting  with  conditions  that  are  compatible  with  the  model  s 


climate  (see  also  Lorenz,  1986).  In  this  run  the  local  planetary 

temperature  T*  varies  vlth  latitude,  but  not  with  longitude  nor  time.  The 

zonal  flow  produced  by  the  cross-latitude  heating  contrast  Is  barocllnlcally 

unstable,  so  that  migratory  waves  develop,  and  transfer  heat  from  one 

latitude  to  another.  At  any  given  latitude,  T  is  generally  not  in 

equilibrium  with  T*.  and  the  diabatic  effects  represented  by  Fig.  2  are 

balanced  In  the  long  run  by  the  advective  temperature  changes. 

The  total  range  of  To  in  Fig.  3  Is  somewhat  greater  than  2.0  K.  Local 

fluctuations  of  this  magnitude  would  be  unspectacular,  but  the  fluttuatlons 

In  Fig.  3  are  of  global  averages,  and  they  seem  to  be  at  least  comparable  in 

magnitude  to  corresponding  fluctuations  in  the  real  atmosphere  during  the 

past  few  centuries  (c/.  Fig.  3  of  Hansen  et  a!.,  1981).  Certainly  an  almost 

monotone  increase  or  decrease  in  the  real  atmosphere  comparable  to  those 

beginning  in  Fig.  3  at  years  32,  114,  or  351  would,  after  progressing  for  a 

number  of  years,  be  interpreted  by  many  investigators  as  a  climatic  change, 

^0 

and  a  cause  would  most  likely  be  sought.  In  M  the  cause  of  these  changes  is 
purely  internal;  the  external  forcing,  represented  by  T*.  undergoes  no 
change  with  time. 

Formally,  the  sensitivity  o(T,T*)  of  T  to  T*  is  given  by  dT/dT\  where 

0*0 

the  differentiation  is  performed  under  conditions  of  constant  F.  We  may 

*0  ^ 0 

also  speak  of  the  sensitivities  o(F.T')  and  o(F,T)  of  the  rate  sof  diabatic 

•*0  ^0 

heating  to  T*  and  T;  they  are  given  by  3F(T’,T)/9T#  and  9F(T’,T)/aT.  It  is 
evident  that  o(F,T)  =  -o(F,T*)/o(T,T’). 

We  now  propose  that  the  long-period  fluctuations  exhibited  in  Fig.  3 
are  related  to  the  high  climatic  sensitivity  o(T.T’)  and  its  restriction  to 
a  limited  range,  which  assures  us  of  a  highly  variable  sensitivity,  although 
the  high  heat  capacity  of  the  oceanic  mixed  layer,  set  at  ten  times  that  of 
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the  atmosphere,  also  plays  a  role.  In  brief,  the  internal  advective 
processes  do  not  alter  To,  but  they  can  alter  the  values  of  T  at  two 
different  locations  by  equal  and  opposite  amounts.  If  some  sort  of 

AT 

equilibrium  has  been  established,  and  if  the  values  of  o(F,T)  at  these 
locations  are  different,  the  alterations  in  the  diabatic  effects  will  not 
cancel,  and  the  equilibrium  will  be  upset.  Since  o(F,T‘>  is  not  highly 
variable,  large  variations  in  the  sensitivity  o(T,T‘)  imply  variations  in 
o(F,T). 

To  support  our  hypothesis  we  shall  construct  a  model,  which  we  shall 
call  M,  and  which  will  be  much  simpler  than  M,  but  will  nevertheless  retain 
the  essence  of  the  processes  believed  to  be  Important  for  the  very-long- 
period  fluctuations.  The  model  M  will,  in  fact,  run  four  or  five  orders  of 
magnitude  faster  than  M,  even  though  M  is  "low  order."  We  shall  then 
demonstrate  that  M  can  produce  a  curve  very  much  like  Fig.  3  when  the 
sensitivity  of  diabatic  heating  to  T  varies  considerably  with  T,  but  not 
when  it  is  nearly  constant. 

2.  The  model 

Our  new  model  M  has  a  single  independent  variable,  time  t,  and  two 
dependent  variables,  the  temperatures  Ti  and  Ts  representing  conditions  in  a 
low-latitude  and  a  high-latitude  belt.  The  planetary  temperatures  Ti*  and 
Tt*  at  these  latitudes  are  prespecified.  The  governing  equations  have  the 
form 


dTi/dt  =  F(Ti\Ti)  +  G  ,  (la) 

dTt/dt  =  F(Ti‘,Ta)  -  G  .  (lb) 


where  F  is  a  simple  function  that  represents  the  diabatic  effects  and  is 
supposed  to  approximate  F,  and  G  is  a  random  variable  that  represents  the 
advectlve  processes  and  does  not  directly  alter  the  mean  temperature  To  = 

(Ti  +  T*)/2. 

Certain  key  points  (T*,T)  in  Fig.  2  have  been  designated  by  the  letters 
A,  B,  C,  0.  E,  P.  Q.  The  point  C  is  the  inflection  point  of  the  zero  line. 

Points  B  and  D  occur  where  the  zero  line  has  a  vertical  tangent,  or 
o(F,T)  =  0,  while  A  and  E,  also  on  the  zero  line,  have  the  same  T‘- 
coordlnate  as  C.  Points  P  and  Q,  on  curves  adjacent  to  the  zero  line,  have 
the  same  T-coordinate  as  C.  The  function  F  contains  several  constants  whose 
values  may  be  assigned  at  will,  so  that  F  may  look  rather  like  or  rather 
unlike  F.  We  shall  consider  the  optimum  values  of  the  constants  to  be  those 
that  produce  the  proper  values  of  T‘c,  T*s  -  T’d,  Te  -  Ta,  and  T#  -  Tp;  here 
a  subscript  refers  to  a  designated  point.  However,  we  shall  not  insist  that 
even  the  optimum  values  make  F  a  truly  close  approximation  to  F. 

To  construct  the  function  F,  we  note  first  that  T*c  and  Tc  are  nearly 
equal  in  Fig.  2,  and  we  neglect  the  difference.  We  note  next  that  even 
though  the  climate  sensitivity  o(T,T*>  is  considerably  lower  for  low  values 
of  T  in  Fig.  2  than  for  high  values,  o(F.T')  Is  also  lower  for  low  values 
of  T;  as  a  consequence,  o(F,T),  whose  variations  we  have  presumed  to  be 
responsible  for  the  long-term  variations  of  To,  is  about  the  same  for  low 
and  high  values  of  T,  and  the  principal  contrast  is  between  the  intermediate 
values  of  T.  near  Tc,  and  the  values  that  are  either  high  or  low.  We  will 
not  alter  this  feature  if  we  do  not  recognize  any  departures  of  Fig.  2  from 
symmetry  with  respect  to  a  180°  rotation  about  C,  and  in  particular  if  we 
make  o(F.T’)  a  constant.  A  simple  function  satisfying  our  requirements  is 


F(T*,T)  =  kU(r*  -  T  +  b2r/(  1  +  |  T|  >1  , 


(2) 


where  k  is  an  inverse  damping  time,  U  is  a  temperature  unit,  t  = 

(T  -  Tc)/U  is  a  dimensionless  temperature,  t*  is  defined  analogously  to  t, 
and  b  is  a  dimensionless  parameter  which  affects  the  extent  of  the  region  of 
infinite  sensitivity  of  T  to  T*.  The  use  of  an  absolute  value  in  Eq.(2) 
introduces  a  discontinuity  in  the  curvature  of  the  zero  line  at  C,  but  this 
does  not  appear  troublesome.  Functions  such  as  arctan  t  or  tanh  t  could 
have  been  used  Instead  of  t/(1  +  It!). 

The  points  A,  B,  D,  and  E  exist  when  b  >  1.  Points  B  and  D  then  occur 
where  dr*/dT  vanishes  along  the  curve  F  =  0,  while  A  and  E  occur  where 
t*  vanishes  on  this  curve.  Since  t  <  0  at  A  and  B  and  t  >  0  at  D  and  E, 
we  find  that  r*»  =  -t'o  =  (b  -  l)2,  while  Tt  =  -ta  =  b2  -  1.  It  follows 
that  (te  -  ta)/(t‘b  -  t'o)  *  (b  +  l)/(b  -  1).  This  allows  us  to 
determine  an  optimum  value  for  b,  after  which  an  optimum  value  for  U 
follows. 

Fig.  4  is  similar  in  format  to  Fig.  2,  and  it  shows  the  field  of 
F(T',T),  which  is  supposed  to  simulate  F(T*,T),  when  the  constants  have 
essentially  their  optimum  values.  The  similarities  of  the  two  figures  and 
also  the  differences  are  apparent. 

Because  the  dlabatic  temperature  changes  are  so  slow,  it  will  be 
feasible  to  integrate  Eqs.  (1)  using  uncentered  forward  differences,  with  a 
time  step  Dt  that  is  supposed  to  be  so  large  that  the  advective  effects 
during  successive  steps  are  more  or  less  Independent.  We  shall  therefore 
let 


G  =  cUR  . 


(3) 
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where  successive  values  of  R  are  chosen  at  random  from  a  uniform 
distribution  with  -1  <  R  <  1,  and  c  is  dimensionless.  It  might  have  been 
preferable  to  let  R  have  a  small  negative  mean  value  Instead  of  a  zero  mean, 
since  the  advectlve  processes  should  tend  to  reduce  the  horizontal 
temperature  contrast  T'  -  (Ti  -  T»)/2,  but  we  shall  not  explicitly  Introduce 
this  refinement,  noting  instead  that  adding  a  constant  to  R  would  be 
equivalent  to  changing  Ti*  and  T2*  by  equal  and  opposite  amounts.  In 
completely  dimensionless  form,  the  equations  of  M  reduce  to 


Dti  =  a[Tt*  -  tt  +  b2ti/(i  +  |  Tt|  )1  +  cR  , 


Dtj  =  a[Ti*  -  ti  +  b*Tt/<l  +  |Ta|)j  -  cR  , 


where  a  =  kDt. 


Although  ti  and  T2  are  supposed  to  be  dynamically  coupled,  it  is 
apparent  from  Eqs.  (4)  that  the  value  of  Tt  has  no  effect  upon  the  behavior 
of  Ti,  and  vice  versa.  Hence  two  experiments,  performed  with  the  same 
values  of  ti*  but  different  values  of  Tt*.  and  with  identical  sequences  of 
random  numbers,  would  produce  identical  time  series  for  ti.  We  do  not  feel 
that  this  is  a  shortcoming.  The  advective  processes,  for  which  the  random 
numbers  are  a  substitute,  not  only  affect  Ti  and  Tt  but  are  affected  by 
ti  and  Tt;  a  pair  of  experiments  with  different  values  of  T2*,  and  hence 
different  behaviors  of  T2,  but  with  identical  sequences  of  random  numbers, 
would  therefore  have  no  meaningful  physical  interpretation,  even  though  such 
experiments  might  prove  useful  for  studying  the  properties  of  M. 


When  Ti  >  0  and  Ta  <  0,  ve  can  rewrite  Eqs.  (4)  in  terms  of  To  and 
t\  We  find  that 

Dto  =  a|To*  -  To  +  b*To/(l  +  2t'  +  T'1  -  To2)),  (5a) 

Dti  =  alT”  -  t'  +  b*(T'  +  t'2  -  To2) 

/(l  +  2t'  +  T'2  -  To2))  +  CR  .  (5b) 

With  b  =  0,  which  would  imply  that  o(T,T*)  is  strictly  constant.  To  would 
soon  become  permanently  indistinguishable  from  To*.  With  b  >  0,  however, 

To  will  never  become  steady  (unless  t*o  =  0),  since  it  will  be  affected  by 
t',  which  will  continue  to  vary  because  of  the  random  advectlve  forcing. 

3.  The  experiments 

In  our  numerical  experiments  with  M  we  shall  choose  Dt  =  7.2  days,  or 
0.02  years,  as  a  reasonable  time  interval  separating  essentially  independent 
barocllnlc  events.  In  the  first  experiment  we  shall  use  essentially  optimum 
values  for  the  constants  that  appear  in  F.  From  Fig.  2  (actually  from  the 
numerical  output  used  in  constructing  Fig.  2)  we  observe  that, 
approximately,  Tc  =  Tc*  «  276  K,  Tb*  -  To*  =  1.4  K,  Te  -  Ta  =  35  K,  and 

T«  -  Tr  =  20  K.  It  follows  that  b  =  1.083,  whence  U  =  100.8  K.  Since  the 

values  of  F  on  the  lines  adjacent  to  the  zero  line  in  Fig.  2  are 
+  20  K/year,  1/k  =  6  months,  so  a  =  0.04.  We  shall  let  c  =  0.01,  thus 

implying  that  the  advectlve  processes  can  raise  or  lower  T'  by  as  much  as 

1.0  K  per  7.2-day  time  step — a  not  unreasonably  high  value  (cf.  Eq.  28  of 
Stone  et  al.,  1982). 
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Anyone  performing  experiments  with  Eqs.  4  will  soon  discover  that  the 
general  behavior  depends  crucially  on  the  values  of  TV  and  Ta*.  Since  the 
model  recognizes  distinct  temperatures  at  only  tvo  points,  ve  cannot  expect 
to  produce  the  needed  contrast  between  the  values  of  o(F,T)  at  these  points 
unless  one  temperature  Is  near  Tc,  while  the  other  Is  decidedly  higher  or 
lower.  This  Is  most  likely  to  happen  if  one  planetary  temperature  Is  near 
Tc,  while  the  other  is  not.  Accordingly,  In  our  experiments  we  shall  let 
TV  as  0.1  and  Ta*  *  -0.01,  implying  that  T»*  *  286  K  and  Ta*  =  276  K. 

We  begin  with  the  essentially  optimum  values  a  =  0.04,  b  =  1.08, 
c  =  0.01,  Tc  -  276  K,  and  U  -  100  K.  Fig.  6  is  constructed  like  Fig.  3,  and 
it  shows  the  variations  of  annual  averages  of  T®  produced  by  If  for  400 
consecutive  years.  The  initial  conditions  were  chosen  arbitrarily,  and  the 
solution  was  run  for  20  years  before  the  start  of  Fig.  6. 

The  qualitative  resemblance  of  the  figures  Is  apparent,  even  though  not 
complete.  As  In  Fig.  3,  there  are  many  nearly  steady  Increases  or  decreases 
of  To  that  resemble  climate  changes,  although  there  are  Instances, 
particularly  in  the  third  and  fourth  centuries,  when  the  trends  fall  to 
persist  as  long  as  ten  years.  Near  the  beginning  and  again  near  the  end  of 
the  first  century  there  are  almost  continual  decreases  of  To  lasting  about 
20  years;  the  latter  Is  followed  by  a  more  broken  increase  for  about  35 
years.  The  resemblance  might  be  made  even  closer  by  effectively  stretching 
the  curves,  perhaps  by  altering  the  values  of  the  constants  so  that  what 
took  place  In  400  years  will  require  600  years. 

Further  evidence  that  the  time  scales  displayed  in  Fig.  5  are  nearly 
right  is  provided  by  Fig.  6.  which  compares  100  consecutive  72-day  or  0.2- 
year  averages,  100  consecutive  1-year  averages,  and  100  consecutive  5-year 
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averages.  Only  the  middle  curve  looks  like  Fig.  3;  72-day  averages  are  much 
too  persistent,  while  5-year  averages  are  not  persistent  enough. 

In  Fig.  7  we  present  100-year  runs  with  ba  =  0.6,  0.8,  1.0,  1.2,  and 
1.4,  while  the  remaining  constants  are  as  before.  As  anticipated,  the 
smallest  value  of  b*  does  not  produce  large-amplitude  fluctuations  of  To. 
and  the  fluctuations  that  do  occur  do  not  seem  to  favor  long  periods.  The 
amplitude  is  greatest  when  b  =  1.0,  where  infinite  sensitivity  of  T  to  T*  is 
Just  appearing.  It  may  seem  surprising  that  the  amplitude  is  smaller  again 
when  b*  =  1.4,  where  infinite  sensitivity  is  well  established,  but  we  can 
explain  this  by  noting  that,  in  the  figure  that  would  replace  Fig.  4  when 
b*  -  1.4,  with  curves  that  are  more  strongly  S-shaped,  the  zero  line  and  the 
vertical  line  T*  =  Ta*  ~  275  K  would  intersect  at  a  value  of  T  far  below  Tc, 
with  a  moderate  value  of  o(T,T*).  There  is  also  a  second  stable 
equilibrium  value  of  Ta  when  ba  *  1.4,  considerably  higher  than  Tc;  our 
experiments  were  performed  with  Ta  oscillating  about  the  lower  equilibrium 
value. 

We  conclude  from  these  experiments  that  contrasting  values  of  o(F,T), 
the  sensitivity  of  the  rate  of  diabatic  heating  to  temperature,  are 
essential  if  long-period  variations  of  To  are  to  be  produced  by  M,  or  by  a 
model  that  Includes  only  those  physical  processes  simulated  by  M. 


4.  Conclusion 

We  have  seen  that  a  low-order  model  of  a  moist  general  circulation  can 
generate  long-term  variations  of  the  globally  averaged  temperature,  bearing 
a  resemblance  to  climatic  fluctuations,  even  when  the  external  influences 
undergo  no  variations  with  time.  The  variations  depend  upon  a  cloud-albedo 
feedback  process  that  develops  In  the  model.  We  have  also  seen  that  a  much 


simpler  model,  in  which  the  barocllnlc  processes  are  represented  by  a 
sequence  of  random  numbers,  can  generate  similar  long-term  variations. 

Our  point  is  not  that  long-period  fluctuations  can  be  produced  by  a 
perhaps  unrealistic  cloud-albedo  feedback  process.  It  Is  that  they  may  be 
produced  by  any  process  or  combination  of  processes  that  can  make  the 
sensitivity  of  diabatic  heating  to  temperature  highly  variable. 

In  effect,  we  have  uncovered  an  apparently  overlooked  mechanism,  to  be 
added  to  a  long  list  of  known  or  suspected  mechanisms,  for  producing  long- 
period  atmospheric  variations.  Whether  this  mechanism  is  present  in  the 
real  atmosphere  or  the  real  climate  system,  and.  if  it  is,  what  the  specific 
processes  that  enter  the  mechanism  may  be,  we  are  not  presently  in  a 
position  to  say. 
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Figure  Captions 


Fig.  1.  The  equilibrium  values  of  the  sea-level  temperature  T.  in 
degrees  K.  corresponding  to  values  of  the  local  planetary  temperature  T‘,  in 
degrees  K,  for  model  M,  when  the  half-life  for  clouds  is  18  hours  (left 
curve).  36  hours  (center  curve),  and  64  hours  (right  curve). 

Fig.  2.  The  dlabatic  heating  rate  F(T‘.T).  in  K/year,  corresponding  to 
the  indicated  values  of  T*  and  T,  in  degrees  K,  for  model  M.  when  the  half- 
life  for  clouds  is  36  hours.  The  points  A.  B,  C.  D,  E.  P,  Q  are  key  points 
used  in  calibrating  the  approximation  F  to  F. 

Fig.  3.  The  variations  with  time  t,  in  years,  of  annual  averages  of 
the  globally  averaged  sea-level  temperature  T*.  in  degrees  K,  produced  by 
model  M.  for  400  successive  years. 

Fig.  4.  The  function  F(T\T),  in  K/year,  serving  as  an  approximation 
in  the  model  M  for  the  function  F  shown  in  Fig.  2. 

Fig.  6.  The  same  as  Fig.  3,  but  for  the  model  M,  with  Tc  =  276  K, 

U  =  100  K,  a  =  0.04,  b  *  1.08,  and  c  =  0.01. 

Fig.  6.  The  variations  with  t.  in  years,  of  0. 2-year  averages  (top 
curve),  1-year  averages  (middle  curve),  and  6-year  averages  (bottom  curve) 
of  To,  in  degrees  K,  produced  by  11,  with  the  conditions  of  Fig.  6. 

Fig.  7.  The  variations  with  t,  in  years,  of  annual  averages  of  T*.  in 
degrees  K,  produced  by  U,  with  the  conditions  of  Fig.  6,  except  that 
bs  *  0.6  (top  curve),  0.8  (next  curve),  1.0  (middle  curve),  1.2  (next 
curve),  and  1.4  (bottom  curve). 
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